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Abstract 

We numerically solve the one dimensional Boltzmann equation of the neutrino and anti-neutrino 
transport in accretion disks and obtain the fully energy dependent and direction dependent neu- 
trino and anti-neutrino emitting spectra, under condition that the distribution of the mass den- 
sity,temperature and chemical components are given. Then, we apply the resulting neutrino and 
anti-neutrino emitting spectra to calculate the corresponding annihilation rate of neutrino pairs 
above the neutrino dominated accretion disk and find that the released energy resulting from the 
annihilation of neutrino pairs can not provide sufficient energy for the most energetic short gamma 
ray bursts whose isotropic luminosity can be as high as 10 52 ergs/s unless the high temperature 
zone where the temperature is beyond 10 MeV can stretch over 200 km in the disk. We also 
compare the resulting luminosity of neutrinos and anti-neutrinos with the results from the two 
commonly used approximate treatment of the neutrino and anti-neutrino luminosity: the Fermi- 
Dirac black body limit and a simplified model of neutrino transport, i.e., the gray body model, and 
find that both of them overestimate the neutrino/anti- neutrino luminosity and their annihilation 
rate greatly. Additionally, as did in Sawyer (2003), we also check the validity of the two stream 
approximation, and find that it is a good approximation to high accuracy. 
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I. INTRODUCTION 



Short Gamma ray bursts (SGRBs) is one of the most energetic phenomena in nature 

[inn 

releasing energy as much as 10 52 ergs in less than 2 seconds [1][2][3]. Until now the central 
engine of SGRBs is still an open question though numerous models have been proposed 
for decades and one of the most commonly accepted model is the neutrino dominate accre- 
tion flows (NDAFs) involving a hyper-accreting black hole with mass accretion rate about 
0.1 ~ 10 M sun /s. The characteristic mass density is about 10 9 ~ 10 12 g/cm 3 and the charac- 
teristic temperature is about 10 10 ~ 10 11 K in the inner part of NDAFs. Therefore, photons 



are almost completely trapped and only neutrinos/anti- neutrinos can escape 
carrying away the gravitational energy of the accreted gas ( 4]j5]j6]j7]0|ji 



rom the disk 
A). Under 



some extreme conditions with high accretion rate (> 1 M sun /s), even the neutrino opacity 
in the disk cannot be neglected, either. When the inner part of the disk is not transparent 
to neutrinos, Di Matteo et al. [6[ investigated the neutrino transport in the disk under two 
main approximations, the first one is that they applied the energy-averaged cross-sections of 
neutrino-matter interaction, which is similar to the Rosseland approximation for photons; 
the second one is the two stream approximation of the angle dependence of the neutrino 



distribution. Sawyer 12] checked the above two approximations by solving the full Boltz- 
mann equation directly, his results have shown that the two stream approximation is good 
enough, but the energy-averaged approximation is only accurate to order. 

Obviously, it is vitally important to solve the neutrino transport problem if we want to ob- 
tain a more precise and more realistic neutrino emitting spectrum, i.e., the energy dependent 
and direction dependent spectrum, which sensitively determines the finial annihilation rate 
of neutrino pairs. Unfortunately, few works focusing on neutrino transport and neutrino 
spectrum in NDAFs have been done and the two common approximate treatment of the 
neutrinos/anti-neutrinos luminosity in the previous works are as follows: first, the neutrino 
emission is integrated over the volume of the disk neglecting the absorption in the neutrino 
transparent cases and assuming the neutrino spectrum to be the Fermi-Dirac black body 
limit in the neutrino opaque cases: f(p) = 1 / (exp(p / kT) + 1) (we denote the assumption 
as Fermi black body spectrum for short in the following sections), where p is the energy of 
neutrinos/anti- neutrinos and T is the local temperature of the disk second, a simplified 
neutrino transport model for both neutrino transparent and opaque cases p, i.e., a gray 



body model for neutrinos/anti-neutrinos spectrum (f(p) = b/(exp(p/kT) + 1), where b is 
a function of neutrino opacity) is applied, please also see Eq.(14) for reference [10 |. Mean- 
while, the neutrino transport problem has been investigated in the atmosphere of neutron 
stars and in the spherically symmetric supernova. In the atmosphere of neutron stars, it is 
supposed that the energy flux and lepton number flux are fixed as boundary conditions and 
are conserved through the atmosphere (Schinder and Shapiro the result shows that 

the only effect of the atmosphere is to soften the neutrino and anti-neutrino spectra because 
the neutrino cross section is roughly proportional to the square of the neutrino energy, and 
therefore the high energy neutrinos are much easier to be absorbed. Neutrino transport 
problem in spherically symmetrical supernovas has also been investigated by Burrows et al. 



151 ] whose result shows that the spectra of neutrinos and anti-neutrinos are neither in the 
form of the Fermi black body spectrum nor the gray body model spectrum. Though the 
atmosphere of neutron stars and supernova are not the same as NDAFs, their results and 



the results of Sawyer [12[ strongly imply that the spectra in the form of the Fermi black 
body spectrum or the gray body model spectrum are not guaranteed in arbitrary situation. 

Although the neutrino transport in accretion disks was investigated by solving the full 
Boltzmann equation directly by Sawyer [12), the main purpose of Sawyer [12] is to check 



the validity of two stream approximation and the energy-averaged approximation, while 



the author did not give the explicit neutrino spectrum. In Ref. 12|, the author also made a 
simplification that the driving chemical potential /i eq = fi c +fi p —fi n (/i c , /i p , fi n is the chemical 
potential of electron, proton and neutron respectively) of neutrinos is much smaller than the 
temperature in disks and can be neglected. However, this simplification will smear the 
difference between neutrinos and anti-neutrinos which is rather important in determining 
the resulting neutrino and anti-neutrino spectra and the final annihilation rate of neutrino 
pairs. 

In this paper, we study the spectra of neutrinos and anti-neutrinos by solving the full one 
dimensional Boltzmann Equation in an infinite and homogeneous disk on the equator plane 
of a central compact object. All the physical quantities vary in the vertical direction and 
are symmetrically distributed about the equator of the disk z = 0. Therefore, the boundary 
conditions for neutrino distribution function /(p, z) can be written as /(p, 0) = /(— p, 0) and 
/(p, H) = for p z < ,where p is the momentum of neutrinos and H is the upper boundary 
of the disk. In §11, the Boltzmann equation for neutrinos and the absorption/emission 



coefficient of neutrinos are presented. Our numerical methods for solving the Boltzmann 
equation are briefly given in §111. In §IV, given the conditions of the disk, we calculate 
the emitting spectra of both neutrinos and anti-neutrinos, which are used to calculate their 

n 

annihilation rate above the disk in §V. As did in Ref. [12[, we also check the validity of two 
stream approximation in §VI. Finally, our conclusions and discussions are summarized in 
SVII. 



II. BOLTZMANN EQUATION 



Following the notation of Sawyer 12|, we define fi = cos(8) for the up moving neutrinos 
and /i = — cos(#) for the down-moving ones where 9 is the angle of neutrino moving direction 
to the vertical direction of disk. We also define the distribution function of the up-moving 
and down-moving neutrinos to be f + (z,p,p,) and f-(z,p,/i) respectively, where z is vertical 
coordinate and p is the energy of neutrinos. For the up-moving neutrinos QQQ, 

their 
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-f+(z,p,fM) + - / dnf-(z,p,n) + f+(z,p,n) 
I Jo 



distribution function is determined by 

^ 0/+(w) = Aa [r{nz )^^ p ) - f + ( z ,p,ri)+\ s 

(1) 

and for the down-moving ones, their distribution function is determined by 

Of—(z p u) 1 r 1 

H t^ 2 — = -A a [f cq {T(z),fi eq ,p) - f_(z,p,fi)]-X a -f^(z,p,/i) + -J dfif-(z,p,{i) + f+(z,p,n) 

(2) 

where A a is the absorption coefficient, A s is the isotropic elastic scattering coefficient, in this 
work, we only include the effects of the isotropic and elastic scattering in the scattering term, 
and here / eq = l/(exp ((p — fi cq )/kT) + 1) for neutrinos, / eq = l/(exp ((p + fi eq )/kT) + 1) 
for anti-neutrinos where /i eq = /i c + p, p — /i n , and /i e ,/U p ,/i n is the chemical potential of 
electron,proton and neutron,respectively. 

Specifically, for the neutrino process v e + n <-> e~ + p, the absorption coefficient for 
neutrinos is shown as follow, 

A. = n n a (1 ~ /c)(1 ~ /p) , (3) 
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where a = a [(l + 3g\)/A]({p+ A)/m e ) 2 y 1 — (m e / (p + A)) 2 (l + l.lp/m n ) , n n is the number 
density of neutron and (1 — / e ) (1 — / P )/(l — f„ q ) is the final state blocking and the stimulated 



absorption correction, here f e , f p is the Fermi distribution function of electron and proton, 
respectively. <r is the characteristic neutrino cross section (er = 1.705 x 10~ 44 cm 2 ), A is the 
mass gap between neutron and proton (A = m n — m p =1.29 MeV), and is the axial-vector 
coupling constant (qa ~ —1.23). 

For the anti-neutrino process T7 e + p -h- e + + n, the resulting absorption coefficient is 
written as, 

\ n . (Wc+XWn) (A) 



where a = cr [(l + 3g A )/4]({p — A)/m c ) 2 y 1 — (m c / {p — A)) 2 (l — 7.1p/m Q ), n p is the number 
density of proton and (1 — / e + ) (1 — / n )/(l — fv e ) is the final state blocking and the stimulated 
absorption correction, here / e + , f n is the Fermi distribution function of positron and neutron, 
respectively. 

For the neutrino/anti-neutrino scattering by neutrons: u c +n — > z/ c +n and V e +n — > V e +n, 
the resulting scattering coefficient is as follow, 

^WilMWjLV ( 5) 
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For the neutrino/anti-neutrino scattering by protons: u e +p — > z/ c +£> and z7 c +j9 — > Z'c+p, 
the resulting scattering coefficient is written as, 

A s = a j- (a sin 4 6 W - 2 sin 2 6 W + I±M) ( JL) 2 , ( 6 ) 
where #w (sin 2 #w ~ 0.23) is the Weinberg angle. 

III. NUMERICAL METHODS 

In order to simplify the Boltzmann equation to be more compact, we define two new 
functions F(z,p,fi) and G(z,p,fi) : 

F = f + + f-, (7) 
G = /+-/-. (8) 

So the Boltzmann equation is transformed to be 



dG{z,p,fi) 
* dz 



V°*S^ = -{X. + X.)G, (9) 
A a (2/ cq - F) + A s (jT 1 rf^F - f) . (10) 



The corresponding boundary conditions f + (z = 0) = f-(z = 0) and f-(z — H) — 
are transformed to be G(z = 0) = and F(z = H) = G(z — H). As the first order 
approximation, we neglect the scattering term in Eq.(10), then Eqs.(9)-(10) is a set of the 
standard differential equations with the boundary conditions. After obtaining the first order 
solution F 1 and G 1 numerically, we use the first order results to calculate the scattering term 
in Eq.(10), after that, the second order solution F 2 and G 2 are obtained in a similar way. 
Repeating this iterative process until a convergent solution is achieved, finally, the neutrino 
spectra F(z,p,/i) and G(z,p,/i) or equivalently f + (z,p,/i) and f-(z,p,/i) are obtained. 

IV. NEUTRINO EMITTING SPECTRA 

For simplicity, we consider a characteristic mass density distribution and temperature in 
the inner part of NDAFs: a distribution of exponent-decreasing mass density, isothermal 
temperature, iso-electron fraction in the vertical direction, that is, 



where p c = 10 11 g/cm 3 , z = 1 km, and electron fraction Y c = (n c — n c +)/n^ = (n e — 
n e +)/(n p + n n ) = 0.4 are adopted in our calculation (n e , n e +, nt,, n p , n n is the number density 
of electron, positron, baryon, proton, and neutron respectively). 

The emitting spectra for neutrino and anti-neutrino (j^) 2 f(T,p,fi) are shown in Fig. 1-2. 
In Fig 1, we plot the emitting spectra of neutrinos and anti- neutrinos at the two fixed angles 
to the vertical direction: 84 degrees (nearly horizontal direction) and degree (precisely 
vertical direction) under two typical temperature T = 6 MeV (left panel) and T = 10 MeV 
(right panel). For comparison, the standard spectrum of the Fermi black body is also shown 
in the figure. To investigate the directional dependence of the emitting spectra of neutrinos 
and anti- neutrinos, given the neutrino energy (p = 2kT), the emitting spectra of neutrinos 
as a function of the emitting direction \x = cos 8 are shown in Fig. 2. The left and right 
panels correspond the results under two typical temperature of the disk (T = 6, 10 MeV). 

As shown in Fig 1 and Fig 2, it is evident that the emitting neutrino spectra are neither the 
isotropic ones, nor the standard Fermi black body ones. It is also notified that the neutrino 
and anti-neutrino emitting spectra are quite different, especially when the temperature is 
not very high comparing to the driving chemical potential /z eq , for instance, the gap between 
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FIG. 1. Neutrino spectrum (^) 2 f(T,p, fi) in two fixed directions at the different temperatures 
(left panel: T = 6MeV, right panel: T = 10 MeV). The neutrino spectrum in horizontal direction 
is shown as the solid line, the anti-neutrino spectrum in the horizontal direction as the dot-dashed 
line, neutrino spectrum in the vertical direction as the dashed line, anti-neutrino spectrum in the 
vertical direction as the dotted line, Fermi black body spectrum as the dot-dot-dashed line, where 
the horizontal direction here means the direction for /i = 0.1 or = 84 degrees to the vertical 
direction and the vertical direction means the direction for fj, = 1.0 or in the vertical direction). 




FIG. 2. Neutrino spectrum (^) 2 f(T,p,fi) versus = cos 9 for a characteristic neutrino energy 
p = 2kT at the different temperatures (left panel: T = 6MeV, right panel: T = 10 MeV). 

the spectra of neutrinos and that of anti-neutrinos for the case of T = 6 MeV is much larger 
than that of T = 10 MeV. Both the anisotropic and asymmetric emission of neutrinos and 
anti-neutrinos will make a difference in the final annihilation rate of neutrino pairs, as we 
will discuss in the next section. 
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FIG. 3. Neutrino number flux fi(-^) 2 f(T,p, ju) versus [i = cos 9 for a characteristic neutrino energy 
p = 2kT at the different temperatures (left panel: T = 6MeV, right panel: T = 10 MeV). 



V. ANNIHILATION RATE OF NEUTRINO PAIRS 



The annihilation of neutrinos and anti-neutrinos into the electron/positron pairs above 



the accretion disk is supposed to be the energy budget of the fireba 



deposition rate by neutrino pairs z/ e , z/ e annihilation is given by 16| 



1 of SGRBs. The energy 



17| 



Q[v c v c ) = \- — ^^^^ T dp T d P '(p+p')(p P r f dnf ^7^(1 -cos e) 2 

4 [m e c z Y(hc) b 3 Jo Jo Ja-k Ja-k 

(12) 

where is the angle between neutrino and anti-neutrino beams, and weak coupling constant 
Ci + C 2 ~ 2.34. 

In this paper, for simplicity, we assume that in the region between r in = 10 km and 
r out = 200 km, the mass density is given by Eq.(ll), Y e = 0.4 and temperature is constant, 
then we can precisely calculate the annihilation rate of neutrino pairs, using the obtained 
spectra in the above section. The results are listed as follows: if the disk temperature 
T = 6 MeV, the total luminosity of neutrino and anti-neutrino are L Ve = 2.0 x 10 53 ergs/s 
and L Vc = 0.7 x 10 53 ergs/s, respectively. The annihilation luminosity of neutrino pairs is 
L Uc u c = 2.3 x 10 50 ergs/s, and the annihilation efficiency rj ~ 0.1%, here the annihilation 
efficiency (77) is defined as 77 = L UcVc /(L Uc + Lp c ). If the disk temperature T = 10 MeV, the 
total luminosity of neutrino and anti-neutrino are L Uc = 7.2 x 10 53 ergs/s and L Vc = 8.8 x 10 53 
ergs/s, respectively. The annihilation luminosity of neutrino pairs is L VqT ; c = 1.5 x 10 52 ergs/s, 
and the annihilation efficiency 77 ~ 1%. 
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As discussed above, the neutrino luminosity and their annihilation rate sensitively depend 
on the disk temperature, the change of the temperature from 6 MeV to 10 MeV leads to the 
two orders magnitude in the final annihilation luminosity! In our opinion, the reasons are 
as follows: first, the higher the temperature, the higher the luminosity of the neutrinos and 
anti-neutrinos. Second, the higher temperature, the harder the spectra of neutrinos and anti- 
neutrinos, which results into the larger annihilation cross section and the higher annihilation 
efficiency Third, if the temperature is high, the effect of the driving chemical potential /j, eq 
which causes the asymmetrical emission of neutrinos and anti-neutrinos is suppressed (see 
Fig.2 and Fig. 3). The nearly symmetrical emission of neutrinos and anti-neutrinos also 
contributes to higher annihilation efficiency. 

To support our above arguments, in Fig.3, we show the neutrino number flux //(^) 2 /(T, p, jj) 
(the energy of neutrino and anti-neutrino is taken to be p = 2kT ) as a function of the 
emitting direction [i = cos 8, under the different disk temperature (T = 6 MeV and T = 10 
MeV). It is obvious that the spectrum of neutrinos and anti- neutrinos nearly coincide at 
high temperature (T = 10 MeV), while the difference between the spectrum of neutrinos 
and anti- neutrinos is huge at low temperature (T = 6 MeV). 

At the end of this section, we will check the validity of the Fermi black body spectrum 
and the gray body model spectrum. In Fig 4, we compare the direction averaged neutrino 
number flux F num = (^) 2 Jo dfif(p 1 ^,)fi of full Boltzmann equation and Fermi black body 
spectrum. It is clearly shown in Fig. 4 that the Fermi black body spectrum overestimates 
neutrino luminosity greatly (about 5 ~ 10 times), so if the Fermi black body spectrum is 
adopted to calculate the annihilation luminosity of neutrino pairs the result will overestimate 
about 25 ~ 100 times in this case. Concretely, the Fermi black body spectrum will lead 
to the following luminosity of neutrino and anti-neutrino by integrating the neutrino and 
anti-neutrino spectra Q]: 

L = ^T 4 S (13) 

where a is Stephan-Boltzmann constant, T and S is the temperature and the area of the 
disk respectively. It is easy to obtain that L(T = 6MeV) = 1.46 x 10 54 ergs/s, and 
L(T = lOMeV) = 1.13 x 10 55 ergs/s under the condition we considered. So the Fermi 
black body spectrum overestimates the total neutrino and anti-neutrino luminosity about 
5.4 times and the corresponding annihilation rate about 29 times when T = 6 MeV, and 
it overestimates the total neutrino and anti-neutrino luminosity about 7 times and their 
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corresponding annihilation rate about 49 times. 

Similarly, the gray body spectrum will lead to the following luminosity of neutrinos and 



anti-neutrinos as [6] 10] 



L = b( 7 -*I*)s = (14) 

\8 J (3/4)(t/2+1/V3+1/3tJ ' 

where r = r a + r s is the sum of scattering optical depth and absorption optical depth, r a is 

the absorption optical depth and numerically, 

r s = 2.7 x 10- 7 Tl lPl0 H (15) 

r a = 4.5 x lO- 7 T? lPw H (16) 

where T\\ is the temperature in unit of 10 11 K, pio is the mass density in unit of 10 10 g/cm 3 
and H is the thickness of disk in unit of cm. Under the condition we consider (p 10 = 10, 
H = 10 5 , and T n = 0.7 and 1.16), it is easy to obtain that r(T = 6MeV) = 0.35, r a (T = 
6MeV) = 0.22, b(T = 6MeV) = 0.59 and r(T = lOMeV) = 0.96, r a {T = lOMeV) = 0.6, 
b(T = lOMeV) = 0.83. So the gray body spectrum overestimates the total neutrino and 
anti-neutrino luminosity about 3.2 times and annihilation rate about 10 times when T=6 
MeV, and it overestimates the total neutrino and anti-neutrino luminosity about 5.8 times 
and their annihilation rate about 34 times when T=10 MeV. 

In addition, it is commonly assumed that the emission of neutrinos and anti-neutrinos 
are isotropic and symmetric in calculating their annihilation rate HQ. But, as clearly 
shown in Fig. 2, 3 and 4, the emission of neutrinos and anti- neutrinos are far from isotropic 
(the intensity in different directions can vary 4 ~ 5 times) and far from symmetric (the ratio 
of the intensity of neutrinos and anti-neutrinos can be 2 ~ 3 times). Therefore, it is no 
doubt that the estimation based on the isotropic and symmetric assumption also contribute 
a non- negligible error in the annihilation rate of neutrino pairs. 

Stricktly speaking, we should include the contribution of the annihilation of the /i neutrino 
pairs and r neutrino pairs to the total annihilation luminosity, but there are two main reasons 
to consider only the contribution from the electron neutrino pairs. The first reason is that, 

nn 

as shown in many previous works [e.g. |6|, llOj, Urea processes (p + e — > n + u e , n + e + — > 
p + u e ) dominate the emission of neutrino pairs in the neutrino-cooled disks. Therefore, it is 
reasonable to ignore the luminosity of fj, neutrino and r neutrino pairs from the disk. The 
second one is that the u e + v e — > e~ + e + channel has the bigger rate than the other two 
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FIG. 4. The direction averaged neutrino/anti-neutrino number flux (F num = (^) 2 Jq <i///(p, //)//) 
at the different temperatures (Ze/t panel: T = 6MeV, right panel: T = 10 MeV). For comparison, 
the corresponding result for the Fermi black body is also shown: solid line for neutrino, dash line 
for anti-neutrino, and dot dash line for Fermi black body. 

flavors chanels, since the process goes through a combination of virtual Z° and virtual 
terms, while the two other species have only the only term. 



VI. TWO STREAM APPROXIMATION 



A simpler way to solve the Boltzmann equation in accretion disk is the two stream 
approximation [19] which replaces the full direction dependent distribution by two streams 
with angle cos(#) = ±l/\/3 to the vertical direction. In Sawyer 12j, the author checked 
the validity of the two stream approximation of the neutrino transport in accretion disk, 
and find that this approximation is good enough. With this simplification, the Boltzmann 
equation is written as 

^ df+ d z ' P) = Uf C(l (T(z)^ cq ,p) - f + (z,p)] + l -\ s \f-(z,p) - f + (z,p)l (17) 

^ df ~dz P) = -^[r(T(z),^p) - f_(z,p)) + \\ s [f-{z,p) - f + (z,p)]. (18) 
With the notation we have defined above, under the two stream approximation, the Boltz- 
mann equation is simplified as follows, 

1 dF{z,p,n) 



VS dz 
1 dG(z,p,fi) 
V3 dz 



— (A a + X 8 )G, 

\ a (2r-F). 



(19) 
(20) 
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FIG. 5. The direction averaged neutrino/anti- neutrino number flux F nura resulting from the full 
Boltzmann equation and the two stream approximation at the different temperatures (left panel: 
T = 6MeV, right panel: T = 10 MeV): solid line and dot dash line are for neutrinos resulting from 
Boltzmann equation and two stream approximation, dot line and dash line are for anti-neutrinos 
resulting from Boltzmann equation and two stream approximation correspondingly. 

The boundary conditions are the same as the original ones. The great advantage of the 
two stream approximation is that the original Boltzmann equation which is an integral- 
differential equation is reduced to a differential equation, which simplifies the original Boltz- 
mann equation greatly. Obviously, the disadvantage of the two stream approximation is 
that we lose the directional information of the emergent spectra of emission. 

In Fig 5, we compare the direction-averaged spectra based on the full Boltzmann equation 
and the two stream approximation. As seen from Fig.5, the two stream approximation is a 
rather good simplification to the Boltzmann equation with high accuracy, if we do not care 
the concrete angular dependence of the distribution function, e.g. we only care the amount 
of the total neutrino number flux or energy flux. 



VII. DISCUSSION AND CONCLUSION 



In the previous works on the neutrino luminosity and the annihilation rate of neutrino 
pairs above NDAFs, the commonly used approximations to determine the luminosity of 
neutrinos and anti-neutrinos are proved to be invalid here. Both the Fermi black body 
spectrum [5] and the gray body spectrum 6] 10] will overestimate the total neutrino and anti- 
neutrino luminosity and their corresponding annihilation rate greatly under the condition 
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we consider. In addition, other common assumptions on calculating the annihilation rate 
of neutrino pairs, such as the emission of neutrinos and anti- neutrinos are isotropic and 



symmetric |5j [18| are proved to be invalid at all, so which will contribute extra errors in the 
estimation of the annihilation rate of neutrino pairs. 

According to the energy spectra of neutrinos and anti-neutrinos obtained from the Boltz- 
mann equation, we find that the emission of neutrinos and anti-neutrinos are generally 
asymmetrical. The asymmetrical emission of neutrino and anti-neutrino is an important 
factor that reduces the annihilation rate of neutrino pairs, and is also responsible for the 
chemical evolution Y e (t) in the disk because the asymmetrical emission of neutrinos and 
anti-neutrinos carries the net lepton number flux. This is an important uncertainty that 
has not been self-consistently dealt with in the previous works, due to the fact that neither 
the Fermi black body spectrum nor the simplified neutrino transport model can properly 
deal with the lepton number flux carried by the stream of neutrinos and anti-neutrinos. For 
instance, the unique property of the atmosphere of neutron stars that the lepton number 
flux and energy flux of the stream of neutrinos and anti-neutrinos are conserved through 
the atmosphere, because of its special distribution of the chemical components jl3| jl^j ]. The 
presence of the atmosphere of NDAFs similar to the atmosphere of neutron stars will soften 
the spectra of the neutrino and anti-neutrino inevitably whatever the incident spectrum is 
and leave the neutrino luminosity unchanged, thus will suppress the annihilation rate of neu- 
trino pairs above the disk. But it is well out of the reach of the Fermi black body spectrum 
or the simplified neutrino transport model. So only through strictly solving the Boltzmann 
equation of neutrino transport can we self-consistently deal with the chemical evolution in 
NDAFs, and then precisely determines the emitting spectra of neutrino and anti-neutrino, 
their luminosity and their corresponding annihilation rate. 

From the strict calculation of the neutrino transport and a coarse assumption about the 
distribution of the temperature, mass density and chemical components in the disk, we draw 
the conclusion that a large(r out = 200 km) area of very high temperature(T > 10 MeV) in 
the neutrino-dominated accretion disk is required if the central engine of SGRBs is indeed 
relevant to the neutrino dominated accretion disk. 

Our assumptions on the mass density, the temperature, and the chemical components in 
the accretion disk are too simple, and we will do more detailed works on this problem by 
considering self-consistently the chemical evolution, the energy balance between the cooling 
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and heating process and the fluid mechanical equilibrium. As shown in this work, the two 
stream approximation is accurate enough as long as we do not care the direction information 
in the emitting spectra of neutrinos, so will be a prime choice for dealing with the problem 
of the chemical evolution and the thermal balance in NDAFs. 
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